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Abstract 

We study the problem of the cosmological constant in the context of the multi- 
verse in Lorentzian spacetime, and show that the cosmological constant will vanish 
in the future. 

This sort of argument was started from Coleman in 1989, and he argued that 
the Euclidean wormholes make the multiverse partition a superposition of various 
values of the cosmological constant A, which has a sharp peak at A = 0. However, 
the implication of the Euclidean analysis to our Lorentzian spacetime is unclear. 
With this motivation, we analyze the quantum state of the multiverse in Lorentzian 
spacetime by the WKB method, and calculate the density matrix of our universe by 
tracing out the other universes. Our result predicts vanishing cosmological constant. 
While Coleman obtained the enhancement at A = through the action itself, in 
our Lorentzian analysis the similar enhancement arises from the front factor of 
e lS in the universe wave function, which is in the next leading order in the WKB 
approximation. 



1 Introduction and Conclusion 



One of the major problems of particle physics and cosmology is the smallness of the 
observed value of the vacuum energy, that is the cosmological constant A. We must 
explain why A is many orders of magnitude smaller than the Planck scale pQ. One of the 
most promising attempts to solve this problem is the one based on the Euclidean wormhole 
effect first proposed by Coleman [2]l*l and studied closely by other authors [1H6] [TJ. In this 
paper, we study this mechanism in the context of Lorentzian spacetime. 

As we review in Section [2J a microscopic wormhole effectively induces the following 
bi-local interactions in the partition function of the Euclidean Einstein-Hilbert action, 

c i3 e- 2S ^ I d 4 xdSV^)^9i3J)Oi{x)0 3 {y), 

id 

where S" w h is the semi-wormhole action and x, y represent the points into which the 
wormhole is inserted(see Fig{T]). 



Figure 1: A sketch of an euclidean universe with a small wormhole insertion. The thin 
tube represents the wormhole. 

Summing over any number of wormholes, it turns out that the system is described by 
the following effective action, 



where we have introduced a pair of operators a and a*, satisfying [a,a^] = 1, which can 
be interpreted as a creation (annihilation) operator of a baby universe. The effect of a 
and a) is to make the cosmological constant dynamical, in the sense that it should be 

1 In [3], Banks also discussed the effect of bi-local interaction. In this paper, we mainly follow Coleman's 
argument. 
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integrated as§ 

^universe = / VgdAexp{- J d 4 x^/g(R - 2A)). 

Coleman analyzed which value of A is favored in the partition function. The result is 
that the integral over the metric can be approximated by a 4-sphere solution and leads 
to a factor as follows, 

■^universe ~ / d A 6 A . 



If one include the effect of other parent universes, which are connected with ours through 
baby universes (see FigJ2J), then the multiverse partition function becomes exp(exp(^)). 
From this argument, he claimed that the wormhole effect could solve the cosmological 
constant problem. 




Figure 2: A sketch of an example of the Euclidean multiverse. Parent universes are 
interacting through baby universes. 



What does this result of the Euclidean gravity imply to Lorentzian spacetime? Naively, 
the 4-sphere solution can be interpreted as a bounce solution. Therefore, the action 
associated with the solution is expected to be the amplitude of a universe tunneling form 
nothing to the size of the 4-sphere (see FigJHJ). However, if we computes the tunneling 
amplitude directly by the WKB method as Vilenkin did [14] , we obtain a factor e~£, 
instead of as Coleman did. Therefore, the physical meaning of the 4-sphere solution is 
unclear, and whether or not the mechanism can work in Lorentzian gravity is doubtful. 

In this paper, we analyze the multiverse in the Lorentzian framework, and see that 
the cosmological constant problem is indeed resolved by the baby universes. First, in 

2 In fact, there are many models in which the wave function of the universe naturally becomes a 
superposition of various possible values of the cosmological constant. See, for example, [8] [9] [10] pj] 
|12j [13] . In our argument, however, all universes in the multiverse should have the same cosmological 
constant, which follows naturally in the baby universe mechanism. 
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a=0 time 



Figure 3: The 4-sphere solution can be interpreted as a foliation of 3-spheres whose radius 



expands from zero to 




Figure 4: A sketch of an example of the multiverse. Parents universes emerges with a 
small size e by a tunneling process. In this example, the initial state has no baby universes 
and the final state has two baby universes. 



Section [21 we review the derivation of the effect of baby universes in the Euclidean path 
integral and obtain its Lorentzian counterpart via a Wick rotation (see FigHJ which is 
the Lorentzian version of FigfS]). 

In Section [3], assuming that there are initially no baby universes and that the par- 
ent universe starts from a small size e via a tunneling process, we calculate the parent 
universe wave function (j>E=o(z)i where z is the sizqj of the universe, by using the WKB 
approximation. The result is 

1 [ z i i 

<f) E =o{z) = r — _ r sin( / k E=0 (z)dz + a), 

'■n/2^z^k E =o(z) Jo 



where E = represents the so-called Hamiltonian constraint, which we will discuss later. 
The factor , 1 behaves like -rju for large z and plays an important role for our 

V fc -E=oO) A 

mechanism. 



3 Strictly speaking, z = a 3 /9 has a dimension of volume. However, for the sake of simplicity, we call it 
; 'size" . 
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From this, we can write down the N-universe wave function as follows, 

where A represents eigenvalues of a + a) and is related to A as A = wh A + A , and 
is the probability amplitude associated with creation of an universcl . The point is that 
the quantum state of N-universe is a superposition of various A, or various cosmo logical 
constants A. Then the state of the multiverse is given by 



-'multi/ 

2V=0 

In Section HJ we compute the density matrix of our universe by tracing out the other 
universes. The result is 

p(z,z')~ / dA e~~\fx\ 2 (f)E=o(z'y(f)E=o(z)exp(^== x\ogz IR ), 
Jo V97T 2 A 

where zm is an infrared cut-off of the size of the universes. The density matrix has a strong 

peak around A = 0, which agrees with Coleman's assertion. However, our mechanism is 

essentially different: While the Euclidean analysis obtains the enhancement through 

the action itself, we obtain the above enhancement factor from the front factor i — of 

e tS , which is in the next leading order in the WKB calculation. 

In Section |5l we explain that what our model predicts to vanish is not the present value 

of the cosmological constant but that in the far future. We also discuss the generality of 

our mechanism. 

2 Effect of Baby Universes 

We first review the Coleman's argument on the effect of baby universes [2] (see also |15j). 
We start from the Euclidean Einstein gravity with a bare cosmological constant Aq, 



J Vgexp(-S E ) = J £>^exp(- J d 4 Xy/g(R - 2A 



))• 



4 In fact, there is also a numerical factor coming from the path integral calculation. We will define /i 
more precisely in Section [3] 
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A wormhole configuration effectively contributes to the partition function as follows, 



CijC 



i, 3 



d A xd 4 y^g(x)^g(y)O i (x)O j (y) exp(-S E ), (2.1) 



where cy are some constants and S w ^ is the action of a semi- wormhole. Summing over 
the number of wormholes amounts to the factor 

If we restrict our attention in the case that 0(x),0(y) are the identity operator 1 , the 
interaction becomes 

exp Qe- 2 ** y d 4 x A v^) VFSOj , (2-2) 

where we have absorbed the numerical constant Cy into 5 w h. By introducing an auxiliary 
variable A, this bi-local interaction can be rewritten as a local interaction as follows, 

J dXexp(-e- s - h X J ^xv^M-y), (2.3) 

where the first term in the exponent additively contributes to the cosmological constant. 
Including wormhole configurations is, thus, equivalent to summation over the cosmological 
constant. 

Alternatively, we can express the wormhole effect by using the following Lagrangian 

S eS = S E + e~ s ^{a) + a) J d 4 x^/gjx). (2.4) 

Here, we have introduced a pair of operators a and satisfying [a, a*] = 1, which can 
be interpreted as a creation/annihilation operator of a baby universal! The last term can 
be regarded as a part of the cosmological constant in each eigenspace of a + a* . Although 
(12. 3p and (12 .4p are equivalent, (12. 4p is more convenient to construct the wave function of 
the universe. 

Finally, we obtain the Lorentzian counterpart by a Wick rotation, 

S = J d A x^-g{x){R - 2A ) - e~ s ^{a ] + a) J d 4 x v / -g(x). (2.5) 



5 To understand this formula, one considers an amplitude between the initial and final state both with 
no baby universe (0| exp(e _5wh (a + af) J d A Xyfg)\£t) . By using the Baker-Campbell-Hausdorff formula, 
it is easy to show that this amplitude recovers Eqn. 
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We use this action to study the cosmological constant problem in the context of Lorentzian 
quantum cosmology. 



3 Quantum State of the Multiverse 



Equation (12.51) implies that the quantum state of the whole multiverse including baby 
universes can be expressed as a tensor product of the parent universes' state and the 
state of the space spanned by a and a) . In this section, we first determine the wave 
function of a parent universe <f>E=a{z) with size z. Then we construct the quantum state 
of the multiverse. 

We will make some assumptions about the multiverse state in order to make the 
physical picture as concrete as possible. However, the conclusion we will obtain does not 
depend on these assumptions, as is discussed in Section [5] 

3.1 Wave Function of a Parent Universe 

We first diagonalize a + a'. For each eigenspace of a + a' with eigenvalue A, the action 
becomes 




where A = ~e 5wh A + A . 

We will consider the homogeneous, isotropic and closed universe: 



ds 2 = -N(t) 2 dt 2 + a{t) 2 dQ 2 , 



(3.1) 



where dfl 2 is the metric on a unit 3-sphere. 



Substituting the metric (13. ip . the action becomes 




In terms of z(t) : 



a(tf 



it can be expressed as 



9 ' 




dt N[z 2 /zN 2 - ((9z) 



/3 - 9A^ 



)]■ 
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The momentum p z conjugate to z is given by p z = —zz/N, and the Lagrangian can be 
written in the canonical form, 

La = p z z - NH A , 

where 



Ua(p z ,z) 



1 



A-^Vz - Ua(z)], where U A (z) 



9 1 / 3 9 



A. 



If we take account of the existence of matter or radiation, the potential becomes 

9 ^ 9 , Cmn.t.t. ^rad. 



U A (z) 



9 ^ C ma tt. 

2 r~ 



(3.2) 



(3.3) 



2z 2 / 3 2 Z 2 4 /3 

where C mat t. and C ra rf. are some constants associated with matter and radiation energies. 
The potential is classified into three types, and we will consider the associated four wave 
functions (see FigJS]). 

(A) (B) (C) 





Figure 5: Three types of potential (A)(B)(C) and the associated four wave func- 
tions(a)(bl)(b2)(c) are sketched. The solid line represents the potential Ua(z), and the 
colored dashed line represents the wave function of the universe <j>E=o{z). Case(A):A > 
and the energy contribution from matter and radiation is large. This corresponds to 
our universe. Case(B):A > and the energy contribution from matter and radiation 
is small. Two wave functions (bl) and (b2) correspond to the two classical solutions. 
Case(C):A < 0. 



To quantize this system via path integral, we take the following metric on the config- 
uration space 



\\S 9l 



[IV I 



tfxy/^g^g&SgWStf'* <x J dt(^(5N) 2 + Na(5af 



(3.4) 
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which is invariant under general coordinate transformation, and this metric leads to the 
volume form of the functional integral 

n a 2 5N5a oc USNSz := [dN][dz]. (3.5) 

Collecting these results, the path integral becomes a quantum mechanical system, 

[dN}[dz\[dp z ] exp(i J dt{p z z - NH A )), (3.6) 

where Ha is defined by (13.21) . 

In the rest of this section, we will determine the quantum state of the universe, as- 
suming that the universe initially has a small size e(see Fi^6]), 




z = e 



Figure 6: The path integral (13. 7\\ is defined as a sum over all histories connecting two 
geometries. 



The amplitude between z = e and z = z is given by the following path integral, 

(z\e-^\ e) = / [dp z ][dz][dN}exp(i I dt (p z z - N{t)H A )). (3.7) 
J Jt=o 

z(0)=e, z(l)=z 

By choosing the gauge such that N(t) is constant T, the path integral of N(t) is reduced 
to the ordinary integral over — oo < T < od_|, 

J dT J [dp z ][dz]exp^i J dt (p z z — T%a) j 

z(0)=e, z{l)=z 



/oo 
d r T{z\e~ iTHK \ e) 
-oo 

Cx (z\5(H A )\e) 

Cx (z\5(H A )( [ dE\<t> E }{<f> E \ )|e), 



6 To be precise, we should integrate only positive T if we fix the time-ordering of the surface T, t=0 and 
Tif—i as in Fig|51 However, we take the integration range as — oo < T < oo to obtain the well-known 
Wheeler-Dewitt equation in the path integral formalism. This procedure corresponds to summing over 
the ordering of the two surfaces too. 
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where C is some constant and we have inserted the complete set {\4>e) } defined b; 

( ( f ) ' E \4 >E ) = 6(E-E'), (3.8a) 
KaI^e) = E\4> e ). (3.8b) 

Therefore, by using 4>e(z) = (z\4>e), the amplitude can be expressed as 

C x <^ =o (e)0£ = o(4 (3.9) 

In other words, the quantum state of the single universe that emerged with size e is^ 

Cx^)IM- (3-10) 

We can calculate 4>e{z) in the canonical quantization formalism. By replacing p z — » 
—id/dz in the Hamiltonian (13.21) . Eqn. fl3.8bp becomes 

V^(^ - U A (z))VZ (j> E {z) = Ecj> E {z). (3.11) 

Note that this leads to the Wheeler-DeWitt equation for E = 0. However, we need to solve 
this equation for general E in order to determine the normalization of the wavefunction. 
After short calculation, it is reduced to the following differential equation 

, d 2 



dz 2 



-40z))Vi<M*) = O, (3.12) 



where 



k 2 E (z) = -2U A (z)-^ 

QA _|_ ^Cmatt. %Crad. 

By applying the WKB method to the function ^/z4>e(z), we can easily solve this equation. 

Before that, we recall that there are four types of universe according to the matter and 
radiation energy and A (see the wave functions (a), (bl), (b2), (c) in Figj5j). In order to 



7 Here, we have assumed that the spectrum of the universe is continuous although it does not seem 
true for the Type (b2) or (c) universe since these states are sort of bound states. The justification of this 
assumption is discussed in | Appendix B| 

8 This analysis is similar to that of [4]. 



9 



write down the solution, we must specify to which type of the universe it belongs. As we 
will see in the next section, the relevant quantity for the density matrix of our universe 
is the integral J dz\(f>* E=0 (e)4>E=o(z)\ 2 (see the exponent in the density matrix (14.11) . where 
\i is defined by ( 13.151) ). For Type (a), (bl) and (b2), the integrals are all divergent. As 
compared to Type (a), however, (j)* E=0 (e) is suppressed for Type (bl), and (pE=o(z) for 
large z is small for Type (b2). Further, the integral for Type (c) is finitqj. Therefore, we 
will concentrate on the universe of Type (a) as the first approximation. 

For Type (a), the whole region of z > is classically allowed and the general solution 
is given by a linear combination of 

1 r z 
4>e{z) = exp(±i / dz'kE(z')), (3.13) 



where the normalization is determined from Eqn. fl3.8ap and the detailed calculation is 
given in Appendix A To determine the solution completely, we must specify the boundary 



condition at z = 0. As the simplest example, we require that 4>e{z — 0) vanishj^j The 
solution satisfying this condition is 



1 r 

<J) E =o(z) = r — Tj= rx sin ( / dzk E=0 (z)). (3.14) 

y/n/2y/Zy/k E =o{z) J 

3.2 Wave Function of the Multiverse 

We can construct the wave function of the whole multiverse from the single universe wave 
function 4>e=o (z)- First, we recall that the quantum system of the multiverse has a pair of 
a and a' , and the eigenvalue A of a + plays the same role as the cosmological constant. 
Parent universes in a fixed-A sector feel the net cosmological constant A = |e~ 5wh A + Ao, 
and evolve according to it. Second, to determine the multiverse state, we must specify 
the initial distribution of A. For example, if there are initially no baby universe as in Fig. 
HI then the initial quantum state can be written as / dX e" A2/4 |A) := We will take 

this case for concreteness although it is not relevant for our argument. 



9 As noted before, it is necessary to modify the model such that E can take continuous values for Type 
(c). The detailed discussion is given in |Appcndix B 



10 The boundary condition might be more complex because the behavior in z < e is highly related with 
the unknown dynamics of the singularity. In fact, our discussion is valid even if we choose other boundary 
conditions, which is discussed in Section [5l 

n It might be helpful to regard a + a* as the position operator \f2x of a harmonic oscillator. 



10 



In order to construct the quantum state of the multiverse |0 mu iti), we need the prob- 
ability amplitude of a universe emerging from nothing to the size e, which we denote as 
fi . Together with the factor in Eqn. fl3.9p or (I3.10p . the weight of each universe /i is given 
by, 

H := no x C x 4>* E =o( e )- ( 3 - 15 ) 

A crucial fact is that /x does not depend on A strongly. This is because 4>* E=Q (e) is a 
smooth function of A as is seen from (I3.14p . and C arising from the path measure should 
have nothing to do with A. 

Collecting these pieces together, we find 



l^multi) = l $ ^>' ( 3 - 16 ) 
N=0 

where |$jv) stands for the N-universe state, whose wave function is given by 

<Z> N (z[ N \--- ,z^) = j d\v N <p E=0 (z[ N) )---<p E=0 {z { ^)® e~T|A). (3.17) 

Here, we note that \l and <j)E=o{zi) depend on A via the net cosmological constant A = 
\e~ s ^X + A . 

4 Density Matrix of Our Universe 

In this section, we study the density matrix of our universe and analyze the preferred 
value of the cosmological constant. We will find that the density matrix has a strong 
peak around A = 0. 

Our universe (denoted by z or z') can be regarded as a subsystem of the multiverse. 
Therefore, we can obtain the density matrix of our universe by tracing out the other 
universes (denoted by z± or zf N \ i — 1 • • • N), 

P{ z i z ) — / „ / ^ Li9 N+1 {z,z 1 ,---,z N ) 9 N+ i{Z,Z 1 ,---,z N j 

N=0 J 

J2jn I dX e ~*\v\ 2 <l>E^W4>E^{z)( f dz"\fi(j) E=0 {z")\A 

N=0 ' J ^ ' 



oc dA e A 2 \n\ 2 <j) E=Q {z')*(t) E=0 (z)ey^[ \ dz"\iJ,(f) E=0 (z")\ 2 ) . (4.1) 
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In the last line, we have changed the integration variable from A to A and the factor e 2 
is understood function of A. 

The integral of the wave function can be evaluated by using the solution f)3.14p for 
Type (a) as 

/ dz\<j) E= o(z)\ 2 = [dz 2 sin 2 ( / k E=0 (z')dz')). (4.2) 

J J TTZk E=0 (Z) J 



Actually, this integral is logarithmically divergent since k E= o(z) ~ v9A is constant for 
large z. If we regulate this divergence by introducing the upper bound Zir of the integra- 
tion, which corresponds to the infrared cut-off of the size of universes, then the integral 
becomes 

Z IR I 

dz\(f) E= o(z)\ 2 ~ —=\ogz IR + --- . (4.3) 
V97r 2 A 

Thus, the density matrix becomes 

f°° _\ 2 \u\ 2 

p ~ / dA |/i| 2 e- 2 E=o (z r )*0 E=o (z) exp(-p= x logz IR ). (4.4) 
Jo V9tt 2 A 



As we have discussed in Section I3.2[ \i does not have a strong A-dependence. Therefore, 
the exponential function in (I4.4p provides a sharp peak at A = 0. This enhancement comes 
from the front factor i — of e lS in the WKB approximation, which is essentially different 
from Coleman's original argument, where the enhancement comes from the action itself. 

5 Discussion 

First, we point out an interesting feature of our model. According to the standard model 

of cosmology, the cosmological constant is measured to be a positive and small non-zero 

value. Our model, on the other hand, seems to predict that the cosmological constant is 

1 

exactly zero because of the enhancement factor e va in the density matrix. However, the A 
appearing in the enhancement factor is defined by the asymptotic value of k E= o(z) — > \/9A 
for z — > 00. Therefore, it is not the presently observed cosmological constant but its value 
in the far future that our model predicts to vanish. Thus, our model does not conflict 
with the current observation. 

Second, we make some comments on the assumptions we made in Section [31 and 
explain that our mechanism works in a quite generic condition. The solution in the 
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Figure 7: A sketch of an example of the multiverse. In this case, the initial state has some 
baby universes. 



classical region is generally a linear combination of 

1 f z 

—= ^=exp(±i / dz'k(z')), (5.1) 

and to be concrete, we have chosen the sine- type solution (see Eqn. fl3.14p ). If we consider 
another situation where the universe keeps expanding [H] , we should use the following 
solution^ 

i r 

exp(-z / dz'k(z')). (5.2) 



z 



However, the enhancement still occurs because it comes from the front factor i — , which 
does not depend on the boundary condition. 

Another generalization is about the initial state of the baby universes. In Section [3J we 
assumed that there are initially no baby universes in the multiverse, which is described by 

A 2 

J dXe ~|A). However, our argument is valid even in the case where there are many baby 
universes at the beginning(see FigJT]). If we consider a general superposition J dX /(A)|A) 
as the initial multiverse state, then the density matrix becomes 

dA f 2 (X = A(A))| /U | 2 S =o(^)*^=oWexp(^= x log* /fl ), (5.3) 

V97r 2 A 

which again has the peak at A = unless the function f(X = A (A)) behaves singularly 
around A = 0. 

From these considerations, it seems that the mechanism we discussed in the paper is 
highly generic, and the problem of the cosmological constant is naturally solved in the 
context of the multiverse in Lorentzian spacetime. 



The sign in the exponent can be understood from the relation p z = —zz/N. 
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Appendix A Normalization of the Wave Function 

In this appendix, we check the wave function (13. 13[) satisfies the normalization (I3.8al) . 



poo 

/ dz^iz^Eiz) = 5{E - E') 
Jo 

Substituting the wave function, the left hand side is 



, dz =exp(±i / dz'(k E '(z') - k E (z'))). (A.l) 

o TTZy , k E (z)k E >(z) J 

Note that the delta function can arise from the integral over the asymptotic region z — > oo. 
For large z, k E (z)k E '(z) ~ 9A and k E > ~ k E ~ ^§{E' - E) ~ n=- z (E' - E), where we 
have used k\ ~ 9A + — + ••• . From these, we can check (j A.10 indeed gives, 

f°° 1 1 

/ d(logz)— =exp(±t — (E' - E) log 2 ) = 8(E - E). 
J 7TV9A V9A 



Appendix B Continuous Spectra of Type (b2) and 

(c) Universe 

In Section [3l we assumed that the quantum state of the universe takes a continuous 
spectrum. In this appendix, by considering a model with a homogeneous scalar field rj(t) 
as illustration, we justify that the spectrum of Type (b2) or (c) can be made continuous, 
and the wave function of the universe can be normalized as in Eqn. fl3.8al) . 

Introduction of a spatially homogeneous scalar field rj(t) with the potential V(rj) gives 
the following contribution to the original Hamiltonian %a, 

m " = + \ V W> 
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where p v is the canonical momenta given by p v = ^ = Here we will take V(rf) = 0. 
To show that the eigenvalue equation of Type (b2) or (c) 

has a solution for any energy E, we write the wave function as e tKv f(z) with some real 
number K. Then, it is easy to understand that (1 B.lj) always has a solution by adjusting 
the constant K appropriately depending on the energy, which we denote by K = K E - 

Furthermore, for Type (c), because the wavefunction f(z) is normalizable due to the 
tunneling suppression(see FigjSJ), we have 

J dzdr}<j>* E ,(z,rj) E <i>{z,7j) = J dzdrj e i( - KE - KE ' )v f E ,(z)f E (z) = some finite value xS(E-E'), 

(B.2) 

where the delta-function is obtained by integrating over 77 since 5(Ke> — Ke) oc S(E' — E). 
Thus, we can always make the solution <f>(z, 77) satisfy the normalization condition (j3.8aj) 
by normalizing f(z) appropriately. 
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